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bX) Abstract 

We have examined the possible construction of a stationary, spherically symmetric and 
^ spatially inhomogeneous wormhole spacetime supported by the phantom energy. The later is 

supposed to be represented by the modified Chaplygin gas equation of state. The solutions 
so obtained satisfy the flare out and the asymptotic flatness conditions. It is also shown that 
CO the averaged null energy condition has to be violated for the existence of the wormhole. 
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O 1 Introduction 

^ One of the most exotic geometries that arise as solutions of Einstein field equations is 

c5 the wormhole. A typical two mouth wormhole connects two arbitrary points of the same 

spacetime or two distinct spacetimes. One observes that any typical Schwarzschild space- 
time contains a singularity at r = making it geodesically incomplete. Ellis [T] first 
observed that the coupling of the geometry of spacetime to a scalar field can produce a 
static, spherically symmetric, geodesically complete and horizonless spacetime and thus 
termed it as a 'drainhole' that could serve as tunnel to traverse particles from one side to 
the other. Later on Morris and Thorne [21 [3] proposed that wormholes could be thought 
of (imaginary) time machines that could render rapid interstellar travel for human be- 
ings. While a black hole possesses single horizon which forbids two way travel (in and 
out) of the black hole but this problem does not arise in the absence of horizon for a 
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wormhole. Unfortunately the existence of wormholes require the violation of the most 
cherished energy conditions of general relativity (null, weak, strong and dominant) which 
are in fact satisfied by any normal matter or energy [1]. In particular, matter violating 
null energy condition is called 'exotic matter' [5j. Later it was proposed that wormholes 
could be constructed with arbitrary small quantities of exotic matter [6l [7] . A commonly 
known form of matter violating these energy conditions is dubbed as 'phantom energy' 
characterized by the equation of state (EoS) p = ujp, where p and p are respectively, the 
pressure and the energy density of the phantom energy, with uj < —1. The existence of 
this matter remains hypothetical but the astrophysical observations of supernovae of type 
la and cosmic microwave background have suggested the presence of phantom energy in 
our observable universe [SI E]- It can exhibit itself as a source that can induce an accel- 
eration in the expansion of the universe. The typical size of a wormhole can be of the 
order of the Planck length but it can be stretched to a larger size if it is supported by 
exotic phantom like matter pUl [TT] . The accretion of phantom energy can increase the 
mass and size of the wormhole and hence guarantee the stability of the wormhole [121 [IS]- 
The astrophysical implications of wormholes are not exactly clear but it is suggested that 
some active galactic nuclei and other galactic objects may be current or former entrances 
to wormholes [H]. It has been predicted that wormholes can also produce gravitational 
lensing events [15]. Since wormholes are horizonless, they can avoid undergoing any pro- 
cess of decay like Hawking evaporation and hence can survive over cosmological times. 
But a wormhole may form a black hole with a certain radial magnetic field (a form of 
magnetic monopole) if it accretes normal matter and consequently loses its structure. 

Earlier, Rahaman et al [16] investigated the evolution of wormhole using an averaged 
null energy condition (ANEC) violating phantom energy and a variable EoS parameter 
Lv{r). We here investigate the same problem using a more general EoS for the pressure 
density namely the modified Chaplygin gas. It is well-known that the wormhole spacetime 
is inhomogeneous and hence requires inhomogeneous distribution of matter. This can be 
made by introducing two different pressures namely the radial and the transverse pressure. 
Our analysis shows that the parameters adopted in the equation of state for phantom 
energy have to be tuned such that the radial pressure becomes negative in all directions 
and for all radial distance. This result turns out to be consistent with Sushkov [17]. 

The paper is organized as follows: In the second section, we have modeled the field 
equations for the wormhole spacetime and proposed the methodology that is adopted in 
the later sections. Next, we have investigated the behavior of energy condition ANEC for 
all the wormhole solutions obtained. Finally, the last section is devoted for the conclusion 
and discussion of our results. 



2 



2 Modeling of system 

We start by assuming the static, stationary, spherically symmetric wormhole spacetime 
specified by (in geometrized units G = 1 = c): 



Here /(r) is the 'gravitational potential function' while b{r) is called 'shape function' of 
the wormhole (see Ref. [18] for the consistent derivation of the above metric). The radial 
coordinate r ranges over [ro, oo) where the minimum value Tq corresponds to the radius 
of the throat of the wormhole. If b{r) = 2m(r), the later being the mass, then Eq. (1) 
represents a 'dark energy star' which may arise from a density fluctuation in the Chaplygin 
gas cosmological background [191 [20]. Note that 6(r = Tq) = corresponds to the spatial 
position of the wormhole throat. We shall, in this paper, assume /(r) = constant for the 
convenience of our calculations. This choice, as a special case, is also physically motivated 
and makes the time traveler to feel zero tidal force near the wormhole [IHIES]- A wormhole 
with small |/'(r)| in the vicinity of the throat is likely to be traversable in the sense of 
having low tidal forces. It also makes the wormhole to be horizon-free. 

We take the inhomogeneous phantom energy which is specifled by the stress energy 
tensor: Too = P, Tu = pr, T22 = T33 = pf Here pr and pt are, respectively, the radial 
and transverse component of the pressure while p is the energy density of the phantom 
energy. It represents a perfect fluid (which is homogeneous and isotropic) if Tu = T22 = 
^33 = Pr = P |21J. Note that in the stellar evolution, the difference Pr — Pt creates a 
surface tension inside star which makes it anisotropic. This feature is generically found 
in more compact stars like neutron and quark stars [22J, contrary to normal stars which 
are majorally supported by radial pressure only against gravity. 

The Einstein field equations {Gap = SirTap) for the metric (1) are 
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The energy conservation equation is obtained from T.^^ = 0, which gives 

, 2 2 

Pr + -Pr Pt = 0. 



(5) 
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This equation can be considered as the hydrodynamic equihbrium equation for the exotic 
phantom energy supporting the wormhole. 

Eq. (2) can be written in the form 

dh 

— = Snr^p. 6 
dr 

Let us choose the modified Chaplygin gas (MCG) EoS for the radial pressure [25] 

Pr{r) = Ap{r) - (7) 

Here A, B and a are constant parameters. The MCG best fits with the 3— year WMAP 
and the SDSS data with the choice of parameters A = —0.085 and a = 1.724 |2S] which 
are improved constraints than the previous ones —0.35 < A < 0.025 ^7]. Recently it is 
shown that the dynamical attractor for the MCG exists at u = —1, hence MCG crosses 
this value from either side u;>— lorco'<— 1, independent to the choice of model 
parameters [28]. Generally, a is constrained in the range [0, 1] but here we are assuming 
it to be a free parameter which can take values outside this narrow range, for instance 
« = — 1 as considered below. This later choice a < makes Eq. (7) a combination of a 
barotropic and a polytropic equation of state. 

Let us take the transverse pressure pt to be linearly proportional to the radial pressure 
Pr as 

Pt = npr, (8) 

where n is a non-zero constant. Thus pr is restricted to satisfy Eq. (7) for a given p while 
Pt is arbitrary in nature due to free parameter n. Using Eq. (8) in (5), we obtain 

, 2 2n 

Pr + -Pr - -yPr = 0, (9) 

which gives 

= Cr^(^-'\ (10) 

Here C is a constant of integration. Since the wormhole is supported by a negative 
pressure inducing exotic phantom energy, it yields p,. < if C < and 1 < < cxo in 

order to obtain finite negative radial pressure. Consequently < if < n < oo. Using 
Eq. (10) in (7), we have 

Ap - - = Cr^^"-'\ (11) 

which can be written as 

ylp"+^ - Cr2("-i)p" -B = 0. (12) 

Note that Eq. (12) is a polynomial equation of degree a + 1 in variable p, which does 
not yield solutions for any arbitrary a. We shall, henceforth, solve Eq. (12) for specific 
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choices like a = — 1, and 1. We shall further employ the following conditions on our 
solutions given below [24]: 

1. The potential function /(r) must be finite for all values of r for the non-existence of 
horizon. In our model, this condition is trivially satisfied since /(r) is taken to be a finite 
constant throughout this paper. 

2. The shape function b{r) must satisfy b'{r = Vo) < 1 at the wormhole throat with 
radius Tq, the so-called flare-out condition. 

3. Further b{r) < r outside the wormhole's throat r > Tq. This condition is a direct 
consequence of the fiare-out condition. 

4. The spacetime must be asymptotically fiat i.e. h{r)/r for |r| oo. 
Now we shall consider the three cases for different choices of parameter a: 
Case-a: If a = 0, then (12) gives 

P=f + (13) 



Using Eq. (13) in (6), we get 



Kn = -J 



IT ^ 2n + 1 



+ Ci. (14) 



Here Ci is a constant of integration. Now ^ — > as |r| ^ oo if = 1 and B = —C. But 
here h{r) = constant and hence gives p = which is an acceptable solution and represents 
vacuum (empty space-time) outside the wormhole throat. This corresponds to vanishing 
pressures i.e. Pr = Pt = 0. This vacuum solution requires C = which in turn leads to 
B = 0. Note that condition (4) can also be met if only B = and n < 0. In figures 1 to 
4, we have plotted the ratio b{r)/r against the parameter r. The Fig. 1 shows that the 
ratio declines as r — > oo, although r is restricted to a certain range. The parameter A 
can assume the value in the range —0.35 < A < 0.025 [29], we choose A = 0.025 for our 
work. 

Further, flare-out condition (2) implies 

b\ro) = < 1, (15) 

which gives an upper limit on the size of throat's radius as 

A \^ , , 
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This requires both A > and C > 0. The throat's radius can be obtained by solving 
6(ro) = To which gives 



ro 



A{2n + 1) 
SttC 



[17) 



This quantity is positive if ^/C > and 2n + l > or A/C < and 2n + l < 0. Similarly, 
condition (3) translates into 



r < 



A{2n + 1) 
SttC 



r > ro- 



(18) 



Note that conditions (2) and (3) are satisfied if Ci — 0. 
Case-b: If a = 1, then (12) gives 

which is quadratic in p and gives two roots of the form 



P± 



C^2(n-i) ^ VCV4("-i) + 4AB 
2A 



(19) 



(20) 



These roots are real-valued if the quantity inside the square root is positive while the 
roots will be repeated if it is zero and complex valued otherwise. We next determine the 
shape function b{r) corresponding to these roots by substituting Eq. (20) in (6) to get 



bUr) 



yl(l + 2n) 



2± 



±[2B7rr\l + — \/AAB + C^r^(n-i)r 



(J2 

5-2n 1 9-6n 4A5r^(i-")^ 



l + 2n 
.4(1 -n)^ 



4(1 -n)' 2' 4(1 -n)' 



C2 



|]/[(n - l)(4ASr^ + CV")]:21) 



Here are two constants of integration whereas is the regularized hyper-geometric 
function. Figures 2 and 3 show the ratios b^{r)/r and b^{r)/r versus r, respectively. Both 
ratios decline for large values of r and approach zero, satisfying the asymptotic fiatness 
condition for specific choice of the parameters. 



Case-c: If a = — 1, then (12) yields 

P = 

Use of this in (6) enables us to write 

6(r) 



C 



A-B 



2(n-l) 



2{n + l){A-B) 



.2(n+l) 



(22) 



(23) 
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In figure 4, the ratio b{r)/r is plotted against r, showing its convergence to zero. Further, 
condition (2) implies 

fe'(ro) = < 1, (24) 

which gives the maximum size of the wormhole's throat 

/A — B\ 2^ 

-o<h^) • (25) 



8nC 

In other words, the throat's radius is given by 



2{n + l){A- B) 
8^ 



1 

2n+l 



(26) 



It requires either A — B>0, C>0 and n > —1 or A — B < and C < 0. This later 
choice of parameters is consistent with the ones that are required for < in Eq. (10). 

As we discussed earher, the relativistic energy conditions are satisfied by ordinary classi- 
cal matter but there are some physical processes where these conditions are violated. For 
example, for a black hole evaporation caused by the emission of Hawking radiation |3Uj . 
The quantized fields in the surrounding of black hole produce massive particles carrying 
positive energy density. Due to energy conservation in the whole process, the negative en- 
ergy density is added to the total energy density of the black hole. Consequently the black 
hole loses mass and its horizon shrinks. The energy conditions are also violated when an 
electromagnetic wave is squeezed resulting in the energy density of the wave to become 
negative, zero and positive at certain wavelengths but the averaged energy density of the 
wave remains positive [23j. One observes that the notion of violation of energy conditions 
is quite ubiquitous at the quantum scale. As the quantum effects allow for a localized 
violation of energy condition, there is a limit to an extent by which these conditions can 
be violated globally. In this connection, the 'averaged null energy condition' (ANEC) is 
specified which states that 



T^pk^'k^dX > 0. (27) 

Here T0/3 is the stress energy tensor. A;" is the future directed null vector, 7 is the null 
geodesic and A is the arc-length parameter. In other words, the integrand must be positive. 
In an orthonormal frame of reference, we have A;" = (1, 1, 0, 0), so that T^^k'^k^ = p + Pr- 
We here adopt the ANEC integral from Visser et al [7] to analyze its violation in our 
model: 

I = f(p + p^.)dV = 2 (p + pr)A7Tr^dr. (28) 



The above integral is called the 'volume integral quantifier' |3l]. It is obvious that the 
above integral becomes negative if p + Pr < 0, the violation of null energy condition 



7 



(see also [33]). The violation of ANEC is the requirement for any phantom matter and 
stability of a wormhole. Now we shall take different p and Pr calculated in each of the 
above cases to evaluate /. Our aim will be to find conditions under which J < 0. We 
shall also consider the case of / ^ which suggests the construction of wormhole with 
arbitrarily small amounts of a phantom energy. 



Case-a Using Eqs. (7) and (13) in (28), we obtain 

/ 



87rC(l+A)^,^2n 



A(l + 2n) 

Note that the above integral gives a finite value ii n < —1/2. Hence we obtain 

^A{l + 2ny ' 



(29) 



SvrC 



(30) 



Moreover, the above integral J<0ifl + A>0ory4> —1. Further, A/C < which 
implies either (1) C < and A > or (2) C > and A < 0. Again the former case (1) is 
consistent with pr < 0. Also / — > 0, if either A —>■ or n —1/2. 

Case-b Using Eqs. (7) and (21) in Eq. (28), we get 



47rr((l + 2A)Cr^'' ± rW^AB + CV4("-i)) 



T[32(n 




Q2 

6n 4A5r4(i~") 



2' 4(1 - n) 



AAB + C^r^i-^-^) 

]/[{2n + l){2n 



(31) 



C 



5)(4A5r^ + CV")] 



The ANEC is violated for particular choice of parameters like B = —3, C = —2 and n = 3. 
Note that we have assumed C2± = for the convenience of our calculations. Under this 
choice of parameters, the two integrals I± will be finite. Also figures (5) and (6) show 
the behavior of the integrals /+ < and /_ < 0, respectively. The plots suggest that the 
two integrals I± tend to infinity for large r, so that an infinite amount of ANEC violating 
matter is necessary to sustain these geometries, which is a problematic issue. However 
this problem can be evaded by considering a matching to an exterior vacuum solution 
which gives a thin-shell wormhole solution [32l |33l [35] . Further, the case of /± — > arises 
if n ^ 1, A ^ -1/2 and + AAB 0. 



Case-c Making use of Eqs. (7) and (22) in (28) yields 

87rC(l + A- By+'^'' 



{A- B){l + 2n) 



(32) 
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The above integral is finite if n < — 1/2. Tfierefore we obtain 

^ l + A~B fA-B' 



2n V 8nC 



(33) 



Furtfier, tlie ANEC is violated / < if either C<OotA — B<0. The wormhole is 
supported by arbitrary small amount of phantom energy if A — B 0. 



3 Conclusion and discussion 



In this paper, we have derived three solutions of wormhole by obtaining different forms 
of b{r). This is carried out by employing the modified Chaplygin gas for the pressure and 
using three specific values of the parameter a. It needs to be mentioned that other values 
of a either don't yield any b{r) or if it does exist than the stability conditions 1 to 4 are not 
verified. Hence we have restricted ourselves to these specific cases as shown in the figures 
as well. The solutions so obtained also satisfy the stability conditions. The pressure and 
the corresponding energy density obtained in each case, violate the null energy condition 
p + Pr < and hence the averaged null energy condition is also violated. These conditions 
need to be violated for the existence of any wormhole solution. 

We performed the analysis by taking the wormhole spacetime to be inhomogeneous and 
anisotropic with non-vanishing transverse pressure. The spacetime needs to be anisotropic 
as it was found that considering an isotropic pressure Pr = Pt = P, for /(r) to be finite, 
one cannot construct asymptotically flat traversable wormhole [10]. In our work, we 
represented the radial pressure by the modified Chaplygin gas and the transverse pressure 
to be linearly proportional to the radial one. The MCG has phantom nature with negative 
pressure. Earlier, Lobo [3l] studied the Chaplygin traversable wormhole and concluded 
that the Chaplygin gas needs to be confined around the wormhole throat neighborhood. 
That work was later extended in [35] using the modified Chaplygin gas and it was deduced 
that modified Chaplygin wormholes may occur naturally and could be traversable. Our 
results are in conformity with their results since our solutions meet the criteria of wormhole 
stability and traversability. 

In a recent paper, Gorini et al [36] have presented an interesting theorem which states 
that in a static spherically symmetric spacetime filled with the phantom Chaplygin gas, 
the scalar curvature becomes singular at some finite value of the radial coordinate r and 
henceforth the spacetime is not asymptotically flat. This result apparently forbids the 
existence of wormholes which are required to be non-singular. The theorem is based on the 
assumptions of homogeneity and isotropy of the spacetime. In case of anisotropy {pt ^ 0), 
the above theorem is not applicable and wormhole spacetime appears naturally. For an 
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isotropic and homogeneous spacetime filled with phantom Chaplygin gas, the asymptotic 
flatness can be achieved by cutting the spacetime at some spatial position r = R and 
glued with a vacuum spacetime, in particular, Schwarzschild exterior spacetime can be 
utihzed |3l]. 
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Figure 1: The ratio b{r)/r is plotted against r with C = 3 and for different values of n = 
—2, —2.5, —3, —3.5, —4, —4.5 which correspond to curves in right to left order. 




Figure 2: The ratio h^{r)/r is plotted against r. The parameters are fixed at i? = 6, C = 7 and 
n = 3. 
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Figure 3: The ratio b-{r)/r is plotted against r. The parameters are fixed at B = —2,C = 7 
and n = 3. 
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Figure 4: The ratio h{r)/r is plotted against r for different values of n = —3, —4, —5, —6, —7, —8. 
The parameters are taken B = 1, A = —5 and C = 2 which correspond to curves in right to left 
order 
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Figure 5: The ANEC integral 1+ is plotted against r. The apparent negative values of 1+ show 
the violation of the ANEC condition. The parameters are chosen as B = —3, C = —2 and n = 3 
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Figure 6: The ANEC integral I_ is plotted against r. The apparent negative values of I_ show 
the violation of the ANEC condition. The choice of parameters is the same as in figure 5. 
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